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Abstract 

We consider chains consisting of several identical subsystems weakly coupled by various types of next neighbor interactions. 
At both ends the chain is coupled to a respective heat bath with different temperature modeled by a Lindblad formalism. The 
temperature gradient introduced by this environment is then treated as an external perturbation. We propose a method to 
calculate the heat current and the local temperature profile of the resulting stationary state as well as the heat conductivity 
in such systems. This method is similar to Kubo techniques used e.g. for electrical transport but extended here to the 
Liouville space. 
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1. Introduction 

As a specific topic of non-equilibrium thermodynam- 
ics, heat conduction has long since been of central inter- 
est. Instead of reaching a complete equilibrium state, 
the composite system under some appropriate pertur- 
bation enters a local equilibrium state - small parts 
of the system approach equilibrium but not the whole 
system. 

Within non-equilibrium statistical mechanics the 
theory of linear reponse, originally developed to ac- 
count for electric conductivity, is a very import method 
to investigate dynamical as well as static properties of 
materials [1,2,3,4]. In this context the famous Kubo- 
formulas [5] have led to a rapid developement in the 
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theoretical understanding of processes induced by an 
external perturbation of the system. However, a direct 
mapping of these ideas on pure thermal transport phe- 
nomena (perturbations due to thermal gradients [6]) 
faces serious problems: Contrary to the case of external 
perturbations by an electric field, thermal perturba- 
tions cannot directly be described by a potential term 
in the hamiltonian of the system. Rather, the thermal 
perturbation is introduced by heat baths with differ- 
ent temperatures coupled to the system, thus calling 
for a more detailed description than is needed for elec- 
tric transport. Nevertheless, those methods are often 
used, eventually because of their immediate success in 
describing non-equilibrium processes [1,7,8,9]. 

Recently, the main focus of considerations on heat 
conduction and Fourier's law has shifted towards small 
(one dimensional) quantum systems [10,11]. Typi- 
cally, these systems are chains of identical subsystems 
weakly coupled by some next neighbor interaction. 
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Based on the Lindblad formalism [12] or on techniques 

of quantum master equations [13], heat baths arc then 
weakly coupled to the chain at both ends. It has been 
found that in such systems the appearance of a normal 
heat conduction depends on the type of the interac- 
tion between these elementary subsystems [14]. Most 
quantum mechanical interaction types show a normal 
heat conduction behavior (constant non-vanishing lo- 
cal temperature gradients), whereas for some special 
coupling the local gradient within the chain vanishes 
(divergence of the conductivity, non-normal scenario) . 

Let us now introduce the model system, which we 
are going to investigate in the following based on the 
full nmncrical integration of Liouvillo-von Neumann 
equation as well as on a perturbation theory in Liouville 
Space. 



2. Model System 

The dynamics of the quantum model for heat trans- 
port is given by the Liouville-von-Neumann equation 
(LvN) for open systems 

|p = /:p. (1) 

Thus we have to consider super operators - here C, - 
acting on operators in Hilbert space, e.g. the density 
operator of the system (see [15,16,17]). The complete 
Liouville operator of the open system under consider- 
ation is given by 

£ = £sys-f£l(Tl)+£2(T2) . (2) 

The first term controls the coherent evolution of the 
quantum system defined by the Hamiltonian H: It is 
defined by its action on the density operator p accord- 
ing to 

/:ByBP = -^[^,/3] , (3) 

like for closed quantum systems. The system H is here 
a chain of N identical subunits with n levels each, cou- 
pled weakly by a next neighbor interaction, thus living 
in a Liouville space of dimension /i?'^ . One could think 
of several concrete model systems, for example spin 
models (n = 2), for which the Hamiltonian would read 



iV-l 

(4) 

The first term is the local part of the Hamiltonian, 
whereas the second defines the interaction between the 
subsystems (a'i^^ denote the Paufi operators of the pth 
spin). Choosing Jx = Jy = Jz we get the Heisenberg 
interaction and for Jz = 0, = Jy an energy transfer 
coupling only (XY model). Furthermore, to avoid any 
bias we will often use a random next neighbor interac- 
tion but without disorder (the same random interac- 
tion between different subsystems). 

The chain is weakly coupled to two heat baths, one 
at each end of the system, given by the super operators 
£i(ri) and Lii^Ti) in (2). This bath coupling could be 
realized by standard Lindblad operators [12,14], well 
known from the theory of open systems in quantum op- 
tics. Another equivalent possibility is to derive a quan- 
tum master equation for the model system leading to 
a special bath coupling like e.g. in [13]. 



3. Current and Local Temperature Profile 

We are interrested in the stationary state of the 
above described system, namely of the LvN equation 
(1) in heat conducting scenario. This stationary state - 
a local equilibrium state - defines a temperture profile 
and a heat current. 

As a measure for the temperature T(/x) of a single 
subsystem we use here the local energy of the respec- 
tive system so that < T < 0.5 in units of the local 
level spacing, c.f. [14,18]. This should be appropriate 
for weakly coupled subsystems within the chain (only a 
very small amount of energy is within the interaction) . 
The operator 

^-"loc — -"loc ~ -"loc ' \P) 

measures the local energy difference between two ad- 
jacent subsystems /i and [i + l {n = 1,2, . . . ,N — 1). 

The energy current operator can be derived 

from a discretized version of the continuity equation 
(formulating a current into and out of the site //, re- 
spectively) 
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i[H, H^^^] = jf'^-i''') - . (6) 

Both operators AH^^^''+^\ ^("."+1), will be used 
now to investigate the stationary state. 

4. Perturbation Theory 

Beside a complete numerical investigation of the 
mentioned stationary state, let us try to get the local 
equilibrium state of the system out of a perturbation 
theory. Since the system of the last section is living 
in Liouville space we have to work in this respective 
space, too. 

Firstly, we consider the unperturbed system as our 
chain with both bath systems at the same temperature 
Ti = T2 = T. Since the two baths have exactly the 
same temperature we expect the whole system to settle 
in a thermal stationary state po: This state should sup- 
port neither a heat current nor temperature gradients 
- it is a global equilibrium state with temperature T. 

The eigenvalues and eigenvectors of the unperturbed 
system are given by the eigen equation 

Co\M = h\pj), j = 0,...,n2^-l, (7) 

where the ket-vectors in Liouville space have been de- 
noted as |. . .). A scalar product of vectors in Liouville 
space can be defined by 

(Pi|p,) = iv{ptft} . (8) 

The (unique) stationary state po is also an eigenvec- 
tor of the system with eigenvalue zero, £.o\po) = 0, 
whereas all other eigenvalues have a negative real part. 
This is due to the fact that asymptotically the system 
should enter the unique equilibrium state |po) regard- 
less of which state the system was at the beginning. No 
other eigenvector is able to contribute to the equilib- 
rium state, i.e. all other eigenvectors must be unstable. 

Since the Liouville operator Co is not hermitian, the 
eigenvectors do not form an orthogonal basis, i.e. 

J2\M{M = G (9) 

is in general not the unit operator in Liouville space 
(c.f. [15]). But with the help of the super operator Q it 
is possible to find a dual basis \p-') = G~^\pj) with the 
property 

J2\pjW\ = i- (10) 

j 



The system will be perturbed now by applying a 
small temperature gradient AT. We are interested in 
the properties of the stationary local equilibrium state 
of the system reached because of this perturbation. The 
Liouville operator of the perturbation is thus given by 

AT AT 
£'(AT) = £i(T+^)+£2(T-^). (11) 

The two environment operators are the same as before 

but now with two different temperatures. 

In the case of a chain of two level systems each of 
the two super operators i = 1, 2 of the bath coupling 

at both ends of the system consists of two transition 
processes (in case of finite temperatures) 

Ci{T) = Wl{T)£l + WjiT)£j , (12) 

with the two rates W^^(T) = (1 - T)Ab and Wj (T) = 
TAb (Ab is the coupling strength of the environment, 
T its temperature) and f/, £j are transition operators. 

With this definition the perturbed super operator 
(11) can be rewritten as 

C'{AT) = jCi{T) + £2(T) + , (13) 

with £ = — f ^ +f Z+fj — ^2 ■ The first two terms are just 
bath systems at temperature T like in the unperturbed 
case, we therefore neglect these terms in the following. 
Like in standard perturbation theory in Hilbert space 
we can calculate the state correction Ap in first or- 
der perturbation theory with respect to the stationary 
state of the unperturbed system po (keeping in mind 
that the eigensystem is not orthogonal) as 

pstat = po + Ap = po 2 f. \P^> ■ 

J = l ^ 

(14) 

For a more detailed derivation of this formula see [19]. 



5. Investigation of the Stationary State 

Since the stationary density operator of the system 
is now given by p = po -f- Ap and since we know that po 
does not give rise to any local temperature difference 
and current, the expectation values of the operators 
defined above are determined only by Ap from (14) . We 
thus find for the local internal temperature gradient 



3 



ATAe 



(p^|g|po) 



Tr{Aii-/;-+^)p,} 



(15) 



and for the local current within the system 

= _ ATAb "1^-' (pi^ IV{j('''''+i)p,} . (16) 



The current as well as the local temperature gradient 
depend linearly on the global temperature difference 
of the bath systems. Under stationary conditions the 
current must be independent of /i, = so 

that (16) can be rewritten as 



J = -k'AT 



(17) 



Eigenstates and eigenvalues entering here the global 
conductivity k' depend only on the mean temperature 
of the unperturbed system, not on AT. Based on this 
K as a global property of the system, including its 
contact properties to the environments, let us call (17) 
"external Fourier's Law" . 

Furthermore, combining (15) and (16), we can define 
a local conductivity within the system 



K 



J 



(18) 



implying also k^'^-'^+i' to be independent of the exter- 
nal gradient AT. 



6. Numerical Investigation of the System 

We compare these result with the complete numeri- 
cal solution of the LvN equation (1), here for a Heisen- 
berg spin chain with four spins. In Fig. 1 we show the 
local conductivity of the two central spin systems /v^'^'^^ 
as a function of the external gradient AT. Indeed, we 
find numerically that /t'^'^' according to (18) does not 
depend on AT (dashed line). The exact numerical so- 
lution of the LvN equation shows a weak AT depen- 
dence (solid line). As expected, linear transport and 
our perturbational theory apply for not too large ex- 
ternal gradients AT only. Note that the full range of 
temperatures < T < oo has been mapped here onto 




cq. (1) — 
cq. (18) 



AT 

Fig. 1. Local conductivity k'^'"') in a Hciscnbcrg spin chain 
of 4 spins as a function of the external perturbation AT; 
the solid line refers to the solution of the full LvN equation, 
the dashed line shows (18). 



system number 

Fig. 2. Tciupcrturc profile of a four spin Heisciibcrg chain, 
as calculated by the full LvN equation. Within the chain 
we find a constant gradient and therefore normal heat con- 
duction. 

the interval [0, 0.5] (AT = 0.3 is thus already a large 
gradient) . 

We address now some remarkable features of very 
small heat conducting model systems, based on the nu- 
merical solution of the complete LvN equation (com- 
pare [14]). 

"Normal" heat conduction (Fourier's Law) is asso- 
ciated with a constant but non-zero local temperature 
gradient and thus a finite conductivity ^('''''+1) = « 
independent of site everywhere in the system (see 
[3,18]). In Fig. 2 we show the temperature profile of 
the Heisenberg spin chain for different coupling types 
within the chain. The data are obtained from the full 
solution of (1). The majority of coupling types, the 
Heisenberg coupling and the random next neighbor in- 
teraction, indeed show the normal behavior in the weak 
coupling limit as demonstrated in Fig. 2 ([14,18]). But 
this "normal" transport type does not always show up: 
A spin-spin interaction consisting of an energy transfer 
coupling only (XY model) leads to a vanishing temper- 
ature gradient off the contact regions (sec Fig. 3). 
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system number jj, 

Fig. 3. Tcmperture profile of a four spin XY model, as 
calculated by the full LvN equation. Within the chain we 
find a vanishing gradient and therefore ballistic transport. 



Tmean [AS] 

Fig. 5. Mean temperature dependence of local conductivity, 
of LvN equation. 




AT(2.3) [AE] 

Fig. 4. Fourier's Law in a four spin Heisenberg chain. For 
different mean temperatures the current is linear in the local 
gradient. Conductivity depends on the mean temperature. 

The vanishing gradient in Fig. 3 implies a divergent 
conductivity within the chain k^'^''^"'"^' . Nevertheless 
the current remains finite because of the resistance at 
the contacts, therefore the global conductivity k de- 
fined in (17) also remains finite for this special cou- 
pling type. Therefore we could state that the "exter- 
nal Fourier's Law" is valid even if Fourier's Law proper 
does not apply. 

For the current, in case of a normal transport behav- 
ior, in dependence of the internal gradient we find a lin- 
ear behavior, like proposed by Fourier's law shown for 
three different mean temperatures of the bath systems 
in Fig. 4. Obviously, the local conductivity depends on 
the mean temperature of the bath systems. Further- 
more, we show the dependence of the conductivity on 
mean temperature (see Fig. 5). 

By comparing these solutions of the full LvN equa- 
tion with the results of the perturbation theory we find 
a very good accordance for not to large gradients, as 
expected. The perturbation theory reflects the vanish- 
ing gradient in energy exchange coupling as well as the 
normal linear gradient for Heisenberg chains. Also the 
conductivities are the same as found in the full solution 



7. Conclusion 

We have considered heat conduction in small quan- 
tum systems built up from identical subsystems weakly 
coupled by some next neighbor interaction. By a per- 
turbation theory similar to that introduced by Kubo 
but extended to the full Liouville space of the system, 
we have been able to derive a quantitative equation 
for the temperature gradient and heat currents within 
such systems. These equations depend only on prop- 
erties of the unperturbed system and linearly on the 
strength of the perturbation AT. 

The most remarkable point of the equation for the 
heat current and the temperature profile is the fact 
that the global temperature difference of the external 
bath systems shows up only as a parameter. This is not 
only a numerical advantage, since a diagonalization for 
different global gradients is no longer necessary, but 
also an interesting physical fact: The heat conductivity 
is independent of the external gradient AT, and the 
"external Fourier's Law" is always fulfilled, even if the 
internal gradient of the system is not constant, as long 
as the perturbation theory applies. 

Our approach does not have the problem of intro- 
ducing a potential term into the Hamiltonian of the 
system, like in standard Kubo formulas for heat con- 
duction. The bath systems, modeled by a Linblad for- 
malism, directly define the perturbation in Liouville 
space. Like in standard perturbation theory in Hilbert 
space, the first order correction to the stationary state 
of the system is expressed in terms of transition matrix 
elements of the perturbation operator and the eigen- 
states and eigenvalues of the unperturbed system. Only 
the non-orthogonality of the eigensystem of the un- 
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perturbed system needs a more careful treatment, for- 
mally the equations are very similar. 

Investigations on further aspects of the derived for- 
malism are in progress. We hope to clarify the ques- 
tion of the different transport behavior (non- vanishing 
and vanishing local gradients) under different coupling 
types. 
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